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Abstract 

The Interacting Instanton Liquid Model is used to explore the role of instanton induced dynamics 
in hadron structure. To support the validity of this model in the chiral regime, the quark mass 
dependencies of several properties are shown to agree with chiral perturbation theory, including the 
density of eigenmodes of the Dirac operator and the masses of the pion and nucleon. A quark mass 
m* = 80 MeV emerging naturally from the model is shown to specify the mass scale above which 
the fermion determinant is suppressed, the zero modes become subdominant, and the density of 
quasi-zero modes become independent of the quark mass. 
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I. INTRODUCTION 



The goal of this work is to gain insight into the mechanism by which the structure 
of hadrons arises from QCD. Although lattice field theory provides a powerful tool for 
solving nonperturbative QCD and is now beginning to successfully calculate experimental 
observables with dynamical quarks in the chiral regime, the physical mechanism by which 
hadron structure arises and the dominant degrees of freedom are not directly evident. Hence, 
in this work, we explore the complementary insight that can be obtained from a model 
that focuses on what we believe to be the dominant degrees of freedom relevant to hadron 
structure. In nonperturbative QCD, the way in which gluons interact with quarks depends 
dramatically on the quark mass. In the limit of large mass, heavy quarks move adiabatically 
in a flux tube potential, whereas in the limit of light mass, a variety of evidence suggests 
that the instanton-induced 't Hooft quark interaction plays an important role and provides 
the mechanism for spontaneous chiral symmety breaking. To understand the structure 
of hadrons containing light up and down quarks, we therefore seek to explore the role of 
instantons and their associated zero modes in hadron structure, and do so in the context of 
the interacting instanton liquid model (IILM). 

In the IILM, the QCD path integral over all gluons is replaced by an effective theory in 
which instantons are the effective degrees of freedom, and the gauge fields of the theory are 
those generated by integrating over the positions, color orientations, and sizes of instantons. 
In the context of this model, we would like to understand the mass range in which instanton 
mediated chiral dynamics is manifested and the extent to which it is described by chiral per- 
turbation theory. In particular, we would like to know the mass scale or scales at which the 
continuum fermion determinant is suppressed, at which zero modes become subdominant, 
and at which the density of quasi zero modes becomes independent of quark mass. 

In the present paper, we set-up the formalism to use the IILM 111 to address these ques- 

fin 

tions. The instanton picture (for recent reviews see [2j, [3|) was originally introduced as a 
model for the QCD vacuum based on semi-classical arguments jj| . It was shown that instan- 
tons lead to spontaneous chiral symmetry breaking by introducing strong non-pertrubative 
correlations between fermionic zero-modes localized around the instanton positions. Spe- 



cific 



eatures of the instanton picture have been observed in a number of lattice studies 
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, 9] and there is evidence that chiral symmetry breaking is correlated with smooth 
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lumps of topological charge, whose profile is consistent with that of singular-gauge instan- 



tons jiol . In addition, instanton-induced correlations in hadrons have been studied in 
a number of phenomenological model calculations, where it was shown that the Instanton 
Liquid Model, (ILM), provides a goo d description of the mass and the electro- weak structure 



of pions, nucleons and hyperons 
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20|,|21|. 



Before trusting the IILM to provide useful insight into the role of instantons and their 
associated zero modes, it is important to verify several essential properties. One key issue 
is to verify the chiral behavior of the spectrum of the Dirac operator. By analyzing the 
dependence on the quark mass of the density of eigenvalues of the Dirac operator, p(X). we 



show that, in the chiral limit, the 



perturbation theory (%pt) result |22|: 



ILM results are consistent with the well-known chiral 



lim lim p(A) = Const + 0(A 2 ) 

A— >0 m q ^0 



(Nf = 2). 



(1) 



In addition, we check that at small but-finite quark masses, the IILM generates mass cor- 
rections to the Dirac spectrum that are consistent with those predicted by xpt. To do so, 
we show that from the structure of the IILM Dirac spectrum computed at different quark 
masses, one can predict the power-law infrared divergence in the quark mass of the scalar 
three-point correlator: 
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This result agrees with the prediction of xpt 221 ]: 
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A second important test is verifying that the IILM effective theory can be used to calculate 
hadron masses and that these masses have the dependence on the quark mass expected from 
Xpt. Since correlation functions in the nucleon and pion channels receive contributions from 
a single instanton whereas in the p and A channels the leading contribution occurs for 
two instantons, the nucleon and pion masses depend on the instanton density to first order 
while the p and A are only affected in second order. Hence, we focus here on the nucleon 
and pion masses. Since we have replaced the functional integral over all gluon fields by a 
sum over instanton fields, a key task is to verify that the effective theory is an adequate 



approximation to the full functional by showing that the correlation function decays at 
large Euclidean time like a pure exponential, as it must in a field theory with a transfer 
matrix, and that spontaneous chiral symmetry breaking occurs as it would in the full field 
theory. In the case of the pion, we demonstrate very clear asymptotic decay as a single 
exponential and use the slope to obtain an accurate measurement of the pion mass at each 
quark mass. These data not only have a quark mass dependence consistent with the form 
expected from chiral perturbation theory, but we show that they yield low energy constants 
consistent with those known from QCD. The nucleon is slightly more problematic, in the 
sense that with the Monte Carlo statistics accessible to us, the combined systematic and 
statistical errors are somewhat larger than in the pion case. We show that these data are 
completely consistent with recent lattice QCD calculations extending well into the chiral 
regime and having substantially smaller errors. Fitting the lattice data with xpt again 
yields sensible low energy constants, so that although our IILM data are not sufficient to 
independently determine these constants, since our data is consistent with the lattice data, 
it is also consistent with these low energy constants. 

Having carefully verified the consistency of the IILM results with lattice results and with 
the behavior expected from xpt, we conclude that the low energy effective theory indeed 
gives a useful approximation to the QCD path integral and hence use this model to explore 
the role of instanton induced dynamics. An essential quantity in the IILM is the overlap 
matrix Ty specifing the probability that a quark hops from instanton i to instanton j, and 
we find that a value of the quark mass m* defined to be equal to an average value of Tij is 
a key parameter in characterizing the quark mass dependence of low energy QCD. In the 
IILM, m* = 80 MeV, corresponding to m n = 600MeV, and characterizes the scale for three 
transitions. It is the mass scale above which the fermion determinant is suppressed, the zero 
modes become subdominant, and the density of quasi-zero modes becomes independent of 
quark mass. The significance of these transitions is discussed in the text. 

The paper is organized as follows. In the next section, we review some of the basic aspects 
of the IILM. In section IHIl we present our analysis of the Dirac spectrum and compare our 
results with chiral perturbation theory predictions. In section IIVI we present the details of 
our calculations of hadron masses and compare our results with lattice data. In section |V] 
we perform chiral extrapolations of the nucleon and pion masses and compare the effective 
parameters of the low-energy chiral perturbation theory of the IILM with those of QCD. In 
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section IVI} we analyze how in the IILM the chiral dynamics is incorporated in the structure 
of the quark propagator and in the fermionic determinant, and we define and compute the 
mass scale m* that sets the boundary of the chiral regime. All results are summarized in 
section IVIIi 

II. THE INTERACTING INSTANTON LIQUID MODEL 

In the IILM, the QCD path-integral over the gauge field configurations is replaced by 
a sum over the configurations of a grand-canonical statistical ensemble of instantons and 
antiinstantons: 

N++N. N f 

Z QCD ~Z ILM = aTIaTj / II ^id{Pi)e- Si ^\[det{ip + im f ). (5) 

N+,N- + ' ' J i i 

Here, dfi, = d£/jd 4 Zidpi is the measure in the space of collective coordinates, color orien- 
tation, position and size, associated with the single instantons. Quantum fluctuations are 
included in Gaussian approximation, through the semi-classical instanton amplitude d(pi). 
In what follows we will work in the thermodynamical limit and adopt the canonical formu- 
lation with N + = N_ = N/2. 

Since instanton-instanton interactions are important to remove large-sized instantons 
from the vacuum Q], the partition function (jSJ) describes an interacting liquid, rather than 
a dilute gas of pseudoparticles. The corresponding interaction action is defined as 

S int = S tot -(N + + N^S , (6) 

where S to t is the total classical action corresponding to a given ensemble configuration and 
So = is the individual single- instanton action in the dilute-gas limit. 

In the IILM, the interaction action is approximated by a pure two-body instanton- 
antiinstanton term, which only depends on the relative coordinates of two pseudo-particles: 

N/2 

Sint — ~ Sintjl i J) (7) 

Such a two-body action is calculated classically, by first determining the total action Stot 
from the gauge potential A^I, A) corresponding to the instanton-antiinstanton pair and 
then subtracting the free contribution 2Sq. Since no exact instanton-antiinstanton pair 
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solution A^IjA) to the classical Yang-Mills equations of motion is known, we rely on the 
stream-line construction, in which the classical action is minimized in all the directions 
except along the collective coordinate describing the separation between the two instantons 



23l | . The analytic instanton-antiinstanton gauge potential in the streamline Ansatz is given 



by 



2 

Xv , o T>abb P 



- 2rj»v o 7 n 2\ + 2R ^~T„2^2 , „2/\\ ( 8 ) 



'x 2 + p 2 X ^ A x 2 (x 2 + p 2 / A) 

where A is the streamline conformal parameter defined by 

„2 \2 



X ~ 2o,o, + \ \~oJo~ 2 l ) (9) 



R = \zi — za 



and 



R ab = -Tr[U^r a Ur b ) (10) 



represents the relative color orientation of the two instantons. 

A major shortcoming of the streamline construction is that, in the most attractive color 
orientation channel, the interaction action smoothly approaches Si nt = —2Sq at short dis- 
tances. This means that the instanton-antiinstanton pair tends to annihilate. The resulting 
gauge field configuration is characterized by a vanishing topological charge and a weak gauge 
field potential, and therefore it corresponds to a perturbative fluctuation. Since perturbative 
contributions cannot be treated consistently in the present approach, Shuryak and Schafer [l| 
suggested removing them by introducing a purely phenomenological short-range repulsive 
core: 

Score = ^4 1^1 ' (H) 

where u is a color orientational factor. Such a term provides a cut-off to the momentum that 
can be exchanged through the instanton field, hence restricting the region of applicability of 
the approach to the non-perturbative sector characterized by momenta of the order p < 1/p, 
where p is the average instanton size. In the language of effective field theory, this repulsive 
core plays the role of a counterterm, parametrizing the ultraviolet physics. However, it 
should be stressed that our hard-core is not derived from a systematic contruction and 
therefore introduces model dependence in the calculation. The coefficient A in (1111) controls 
the strength of the repulsion and is the only phenomenological parameter of the model. In 
our calculations we adopted the value A = 128 suggested by Schafer and Shuryak [l|. 



Once the interaction action is defined, it is possible to compute the main properties of 
the ensemble, the instanton density n, and instanton size distribution d(p), by minimiz- 
ing the ensemble's free energy numerically. One can then calculate arbitrary Euclidean 
Green functions by performing Monte Carlo averages over instanton configurations, using 
the Metropolis algorithm. As in analogous lattice QCD calculations, one must set the scale 
in physical units by matching one dimensionful quantity. 

In order to compute correlation functions involving quark field operators, we follow the 
same prescription adopted in lattice simulations, i.e. we first explicilty integrate-out the 
fermion fields and then compute Monte Carlo averages of the resulting Wick contractions. 
Such averages are performed using configurations which are obtained by an accept/reject 
Metropolis algorithm in which the contribution of the fermionic determinant is included in 
the Bolzmann weight, corresponding to unquenched simulations. 

Unlike in lattice QCD simulations, where the fermionic determinant and the quark prop- 
agators are evaluated using a purely numerical algorithm, in the IILM these quantities have 
a semi-analytic representation. This property allows us to identify the physical content of 
each term and establish a connection with the corresponding non-perturbative quark dy- 
namics 2J. The fermionic determinant in a given instanton gauge field background is 
factorized into a contribution of near zero-modes and a contribution arising from non-zero 
modes: 

Det A (# + m q ) = Bet Azm x Det' Anzm . (12) 

The contribution of near-zero modes can be determined exactly, by expanding the Dirac 
operator on the basis of zero-modes of individual instantons: 

Det A2m = Det(T + m g ), (13) 

where Tjj = J d^zip^ \z)iD ^ ^^{z) is the overlap matrix which represents the probability 
amplitude for quarks to "hop" from the instanton / to the instanton J. The non-zero mode 
part of the fermionic determinant is approximated with the product of the non-zero mode 
contributions of each individual instanton: 

N++N- 

Det^ m = n (l.Um qPi ) (14) 
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Simlarly, the quark propagator consists of a zero- mode part and a non-zero mode part: 



S(x,y) A = S zm (x,y) A + S nzm (x,y) A . 



(15) 



The zero-mode part of the propagator S zm (x, y) A dominates in the low-energy regime and 
encodes the information on the physics of chiral symmety breaking and axial anomaly satu- 
ration. It is constructed by representing the Dirac operator in the basis of zero-mode wave 
functions ip°(x) of the individual instantons and reads: 



The non-zero-mode part is approximated as the sum of the non-zero-mode propagators in 
the field of the invidual instantons Q]. 

III. THE SPECTRUM OF THE DIRAC OPERATOR AND THE CHIRAL 
REGIME IN THE IILM 

In the previous section, we have seen that the main approximation of the IILM con- 
sists in replacing the QCD path integral with a statistical sum over a classical ensemble of 
pseudoparticles. In this section, we address the question of whether such an approximation 
provides a realistic description of the non-perturbative quark-gluon dynamics associated 
with the spontaneous breaking of chiral symmetry in QCD. To this end we focus on the low- 
virtuality sector of the spectrum of the Dirac Operator Ljt, which encodes information about 
the dynamics associated with the spontaneous breaking of chiral symmetry. Such a connec- 
tion is made explicit in the Banks-Casher relation which relates the density of eigenvalues 
with small virtuality to the quark condensate: 



We recall that this relation holds in a infinite- volume system and in the chiral limit. Our 
goal is to explore finite-mass corrections to the Banks-Casher relation in the instanton model 
and check if such corrections are consistent with predictions derived from xpt. 

In order to set the framework of our analysis, it is useful to briefly review how the Banks- 
Casher relation (flTl) was obtained. The quark propagator S(x,y) A in the fixed background 




(16) 



7Tp(X = 0) 



(17) 
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is expanded in terms of the eigenvalues A and eigenvectors of the Dirac operator as 

S(x,y) A =y — . 18 

^ m q - i\ n 

n H 

Chiral symmetry implies that for every non-zero eigenvalue A with eigenvector ip\ there is 
another eigenvalue —A with eigenvector jsip\. Hence, setting x = y and using the orthonor- 
mality of the basis of eigenfunctions of ip\ we obtain 

*/d^ TV [*(*,,)]. ^ (19) 
Introducing the spectral density p(A) = {^2 n S(\ — A n )) we can rewrite (fT9l) as 

(<l<l} = -2m q r dX -T^V2 (2°) 

J0 m q + * 

At this point we take the limit V — > oo and m — > 0. We notice that the order in which 
the two limits are taken is crucial since we only have a finite quark condensate in the 
thermodynamic limit: 

(qq) = lim lim -2m„ [ dX P } X \ = -vrp(A = 0). (21) 

Eq. (ITT]) can be interpreted as the lowest-order term in the Taylor expansion of p(A) near 
the origin, in the chiral limit. The next order in A was derived by Smilga and Stern 22] 
using xpt: 

1 B 2 (N 2 t - 4) 

where Nf is the number of flavors, B and fo are the constants that appear in the lowest 
order xpt Lagrangian 

C = / 2 (i Tr[d^Ud^) + IB Tr[MU ] + UM*]), 

U = exp (i—fo), and M is the diagonal quark mass matrix. 

We stress the fact that the result (|22j) is valid only in the limit in which the quark mass 
is much smaller than the typical eigenvalue saturating the spectral integral ( !20l) . m q <C A. If 
the quark mass is of the order of the typical eigenvalue A, then mass-dependent corrections 
to the spectral density profile are expected to appear. Hence, the question arises if there 
exists a range of masses for which the xpt predictions are accurate, yet mass corrections to 
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FIG. 1: Spectral density of the Dirac operator for N c = 3 and Nf = 2. The lower panel shows 
an enlargement of the upper panel in the region < A < 0.3 and the solid lines represent best fits 
using p(X) = c\ + C2 x m^A a_1 . 

( 122]) become important. In order to address this question, Smilga and Stern considered the 
3-point scalar correlator: 

K abc = J d 4 xd A yd 4 z(0\S a (x)S b (y)S c (z)\0), (23) 
10 



where S a (x) is the scalar density operator, S a (x) = q(x)r a q(x). By proceeding in the same 
way as for the quark condensate they obtained an expression for K abc in terms of the density 
of eigenvalues, p(A): 

Rabc = _ d a b c / (A V I J d \. 24 

Jo {m 2 q + X 2 Y 

In the limit in which all mass corrections to (I22p can be neglected, one can obtain a prediction 
for the 3-point scalar correlator by substituting the expansion of p(A) given in (|22|) in (|24|) 
and performing the integral. The constant term in p(A) does not contribute, while the linear 
term gives rise to a power-law singularity in the pion mass: 

B 3 (N 2 f -12) 1 

K ahc = -±J——±— oc l/{ml). (25) 
On the other hand, evaluating the same correlator to lowest-order in xpt one finds: 

*~ = Sp?^' (26) 

Note that, although in both cases one finds the same power-law dependence on the pion 
mass, the numerical coefficient turns out to be different 

jV|-4^(iVf-12)/2 (27) 

The explanation of the mismatch is that the characteristic eigenvalues, A, saturating the 
power divergent integral (T24T) are of order m g , so we are out of the range of validity of (f2~2l) . 
In order to reproduce the correct xpt prediction for the scalar 3-point correlator such mass 
corrections have to be included. 

It is particularly interesting to consider the case Nf = 2, for which the linear term in ( f22l) 
vanishes and the distribution of eigenvalues near the origin is flat. In this case, the power-law 
divergence in the pion mass predicted by xpt must arise entirely from mass corrections. This 
observation can be used to check the consistency of the IILM with xpt. We have performed 
simulations of the Dirac spectrum at different values of the quark mass, ranging from to 
400 MeV. The contribution of the near zero-mode zone to the Nf = 2 spectral density of 
the positive eigenvalues is presented in FigO 

Some comments on these results are in order. When the quark mass approaches zero, 
our results become consistent with the flat trend predicted by fl22|) : 

lim lim p(A) = Const + 0(A 2 ). (28) 
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The value of the constant depends on the normalization of the spectrum and can be fixed 

lis case has been already considered in previous studies 



by imposing the condition ([T 
of Shuryak and Verbaarschot |24j . l25j |. On the other hand, as the quark mass increases, 
the structure of the spectrum changes appreciably. The flatness near A = disappears and 
a peak near the origin develops. In Fig{T] we see that for small A, all curves at different 
masses are very well fitted with functions of the form p(A) — p(0) oc m a \ 1 ~ a . If this form is 
introduced in the spectral representation of the scalar three-point function (1251) . the integral 
can be carried-out analitically leading to results proportional to the inverse quark mass. 

To directly compare the xpt formula for K abc with the results of IILM calculations, it 
is convenient to consider the combination m a Kabc , m . In fact, using the Gell-Mann Oaks 
Renner and Banks-Casher relations, we find the simple xpt prediction 

K abc B{N 2 f - 12) 
m a lim — = —k = const. (29) 



On the other hand, from the spectral representation (pHl) we find 

K abc a f Pmq (\) K-3A 2 ) jabe 
m n — = —mz / dA - — I n^rd , 



Vm,=o(0) q J ~' Pmq=0 (0)(mt + \^- 

which can be computed directly from our IILM points. 

The results are plotted in Figj2j, and we observe that the expected constant behavior with 
m is obtained for m q < 80 MeV. Thus, in this case, the IILM reproduces the structure of 
the Dirac spectrum dictated by chiral symmetry. 

IV. PION AND NUCLEON MASS IN THE IILM 

In this section we present our IILM calculation of the nucleon and pion masses performed 
at different values of the quark masses and we compare with the available lattice results. 
Exploratory estimates of the nucleon and pion masses in the IILM were performed at a 



relatively large quark mass in 26|]. However, such calculations were based on the analysis of 
short-range point-to-point Euclidean correlators and relied on a specific assumption for the 
spectral density (i.e. pole-plus-continuum parametrization) . In the present work, we have 
chosen to avoid any a priori assumption on the spectral density. Instead, we have performed 
the effective-mass plot analysis, which is routinely used to extract the lowest-lying hadron 
masses from lattice QCD simulations. 
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FIG. 2: m„ g °L at different quark mass values in the IILM. 
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To compute the mass of the pion, we have evaluated the correlation function: 



GV(r)= / d 3 x(0|T[j 5 (x,T)j 5 (0,0)|0>, 



where js{x) is the pseudoscalar density operator 



(30) 



f b {x) = q{x)^r a q{x). 



(31) 



To compute the nucleon mass we have evaluated the correlation function 



G n (t) 



d 3 x(0|T[j 7V (x,r)j JV (O,0)P + |0), 



where P+ 



1-74 



is the positive-parity projector and 



(32) 



j a N (x) = e abc u a {x)C lb d b {x)u c {x) 



(33) 



The integration over the final position x in (1301 and fl32l) ensures projection onto zero- 
momentum states. The mass of these hadrons can then be extracted from the plateau in 
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the large Euclidean time limit of the effective mass, i.e. using 



M n/N = lim M e J/ N (r) (34) 
M eff (r) - 1 l n G -/n( t ) ( or.\ 

The point-to-point correlation functions in ( l30i) and ( 1321) were computed by performing 
Monte Carlo averages over instanton ensemble configurations, as described in section [III 
The momentum-projection integrals in (I3"U1) was carried-out using an adaptive Monte Carlo 
routine (VEGAS) with two iterations of 6000 points. For the momentum-projection of the 
nucleon 2-point function (|32|) we combined the results of VEGAS with those obtained using 
a cubic grid with grid spacing 0.1 fm. Statistical errors over ensemble averages were obtained 
using the jackknife technique, with a bin size of 10 configurations. We have used five sets of 
250 independent configurations corresponding to quark masses ranging from 20 to 90 MeV. 
In order to reduce finite volume artifacts, we performed our computations in two boxes of 
size 3.45 3 x 5.9 fm 4 — for the two lightest quark masses — and 2.96 3 x 5.9 fm 4 — for the other 
masses — . We have checked that with such a choice the condition m^L > 5 was always 
satisfied. 

All simulations where performed using the lTFlop/s ECT* cluster, for a total computa- 
tional time of 30000 cpu hours. In Fig. [3] we present pion and nucleon effective mass plots 
for the different values of quark mass. We observe that the all pion effective mass plots 
display a very clean plateau, from which it is possible to unambiguously read-off the pion 
mass from a correlated chi-square fit. On the other hand, the interpolation of the effective 
mass plots of the nucleon is somewhat more problematic. In fact, not only statistical errors 
are larger, but also some of the plots are consistent with a slight slope, in the large Euclidean 
time regime. Hence, one necessarily needs to account for the systematic errors, which can 
be estimated from the smallest and largest mass compatible with the effective mass plot, in 
the nearly flat region r > 1 fm. 

We recall that in the instanton model, the average instanton size p is the parameter 
determining the ultraviolet cut-off and playing the role of the lattice spacing a in lattice 
gauge theory. Hence, we have expressed all dimensional quantities in units of p and then 
determined it from a matching condition. For simplicity, we have chosen the condition that 
the IILM prediction for the nucleon mass should match the CP-PACS lattice result at the 
pion mass m n = 525 MeV. The values of the nucleon masses for m n < 500 MeV are therefore 
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FIG. 3: Pion and nucleon effective mass plots in the IILM at different quark masses 
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FIG. 4: Comparison between nucleon masses obtained with the IILM for three different choices 
of the average instanton size p and those calculated on the Lattice by the CP-PACS [27]] and 
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predictions of the model. The choice of the matching point was motivated by the fact that 
the pion mass ~ 500 MeV represents an upper bound for the mass regime of QCD where 
the dynamics of chiral symmetry breaking is expected to play an important role. 

In Fig. H] we present the results of our calculations for the nucleon mass as a function of 
the pion mass squared for different choices of the cut-off p. Comparing these IILM data at 



different p with the lattice data collaborations 



27 



clearly shows that agreement between 



IILM and lattice QCD predictions is obtained only for a unique value of the cut-off scale, 
corresponding to an average instanton size of 0.32 fm, in excellent agreement with the early 



phenomenological estimate p ~ 1/3 fm [4j and recent lattice results jf], Q, 30, 31, 32]. The 
complete set of values of masses extracted from our IILM calculations is presented in Table [B 
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TABLE I: Pion and nucleoli masses fitted by effective mass plot with x 2 / n df < 1- The quark 



masses are determined at a scale 2 GeV as in 



34] 



m q [GeV] Pion [GeV] Nucleon [GeV] 



0.021 
0.03 
0.05 
0.07 
0.09 



0.300 ±0.004 1.11 



-0.05 
-0.05 



0.360 ±0.004 1.15t{j:o7 



0.460 ± 0.004 1.20 
0.530 ±0.004 1.28; 
0.600 ± 0.004 1.35 



-0.02 
-0.02 

-0.01 
-0.01 

-0.02 
-0.1 



V. CHIRAL EFFECTIVE PARAMETERS IN THE IILM 

In the previous sections, we have shown that the IILM provides a realistic description 
of the microscopic dynamics responsible for chiral symmetry and that it contains pions as 
low-energy vacuum excitations. Thus, we conclude that the effective theory approximates 
the path integral sufficiently well that the IILM can be described by chiral perturbation 
theory, and we therefore ask the next question of how similar the low energy constants are 
to those arising in QCD. Hence, we determine the quark condensate and the pion decay 
constant from the dependence of the pion mass on the quark mass. Note that the numerical 
value of the quark mass in QCD depends on the renormalization scale. In the IILM we do 
not have this freedom since the ultraviolet cut-off scale is provided by the inverse instanton 
size 1/p — 600 MeV. Chiral perturbation theory to 0{p A ) predicts a dependence of the 
form: 

m * = 2m * B °( 1+ 3#b^))- < 36 > 

The chiral scale A is set at the p vector meson mass, which in this model is found to be 



independent on the quark mass, with the value M p = 1 GeV 35j. Then, using the value for 
= 340 GeV" 1 which was extracted from the analysis of the Dirac spectrum ( 1291) . one can 
extract B$. 

The results of the chiral fit to the numerical IILM calculations are shown in Fig. [5j The 
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FIG. 5: Chiral extrapolation of the pion mass, obtained in the IILM using 0(p 4 ) xpt- 



corresponding low-energy chiral coefficients calculated in our model are 



fo = 0.085 ± 0.003 GeV, 
B = 2.43 ±0.02 GeV, 



(37) 
(38) 



corresponding to a chiral condensate of 



(qq) « -(0.259 GeV) J 



(39) 



The fact that these quantities are rather close to the corresponding values extracted in QCD 
implies that the low-energy effective theory of the IILM is indeed not far from that of QCD. 
Note that a recent analysis of the quark mass dependence of the chiral condensate in the 



instanton vacuum leading to comparable results can be found in 



36]. 



A similar analysis can be carried out in the nucleon sector, where several approaches to 



chiral extrapolation have been proposed, for example see 



37 



38 



39 



40 



therein. One can, eg., follow Procura et al. 37], [38( and consider covariant baryon chiral 



41] and references 



perturbation theory (B%PT) with infrared regularization at order 0(p 4 



. In this scheme 
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the nucleoli mass expansion in powers of the pion mass gives the following result: 

M N = M - 4 Cl ml - -f^ml 

3 (9a C2\ 3 f g\ Q „ , „ , A „\,_m 



^ (A) - 64^ { W - 2 ) ~ lilt ~ 8Cl + C2 + bg A J 



4 



256tt/2M 2 
where e£(A) is the finite part of 



^"X (40) 



ei = e'i(A) + ^ (|| - 8 Cl + c 2 + 4c 3 ) , (41) 

and encodes information on the unresolved ultra-violet physics. 

We have performed a fit of the MILC and CP-PACS data in the region < 0.6 GeV, 
supplemented by the nucleon mass at the physical point. To reduce the number of free 
parameters, we fixed f n and qa to their experimental values and c 2 = 3.2 GeV -1 , in ac- 
cordance with the analysis of the ttN scattering in [43j]. The result of our fit is shown 
in Fig. [6J where it is compared with our predictions for the nucleon mass in the range 
300 MeV < Mn < 600 MeV. Although the present statistical and systematic errors in the 
calculated proton mass do not allow us to attempt a direct quantitative determination of 
the effective parameters of B%pt for the IILM alone, the overall agreement between IILM 
and lattice data indicates that they are consistent with those extracted from QCD at the 
level of one standard deviation. 



VI. INSTANTON-INDUCED DYNAMICS IN THE CHIRAL REGIME 

Having established that for < 500 MeV the IILM agrees at a qualitative and even 
quantitative level with several QCD predictions obtained from xpt and lattice calculations, 
we can now use the physical insight from that model concerning the instanton-induced 
correlations to explore the dynamical mechanisms involved in the transition into the chiral 
regime. 

From ffl6|) . it follows that the non-perturbative dynamics associated with the instanton- 
induced near-zero modes becomes parametrically small for quark masses much larger than 
the typical value of the overlap matrix element Tjj. In FigJ7|we have plotted the distributions 
of the maximum overlap matrix elements Tja obtained at different values of the quark mass. 
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1.6 




FIG. 6: Nucleon mass chiral extrapolation 
mined according to the method outlined in 



at 0(p 4 ) in relativistic Bxpt. The error band is deter- 



. The IILM, CP-PACS, and MILC data are denoted 



by the squares, triangles, and diamonds respectively, and the star represents the physical point. 

We note that such distributions are peaked around the value m* ~ 80 MeV. The value m* 
plays a central role in specifying the scale for instanton induced chiral symmetry breaking in 
the IILM. Physically, we expect that for m q 3> m*, corresponding to m n 3> 500 MeV, the 
interactions associated with chiral dynamics become sub-leading. We note that this number 
is consistent with the phenomenological estimate derived in [44|, using the single instanton 
approximation. To further illustrate this transition, in Fig. [8] we show the point-to-point 
pseudoscalar correlation function evaluated in the IILM in the zero-mode approximation 
(zma) in which 



S(x,y) = S free (x,y) + S zm (x,y). 



(42) 



We normalize this full correlation function to the corresponding correlation computed in 
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FIG. 7: Distribution of the maximum fermionic overlap matrix elements in the IILM for different 
quark masses. 

the free massless theory, 

n(r,m ? ) ^ _ (0|j 5 (r)jj(0)|0> IILM—free+zm //fo\ 

n (r) = 4VV^ 

j B (x) = u(x)>y 5 d(x), (44) 

evaluated at a typical non-perturbative scale r = 1 fm for different quark masses, and com- 
pare it to the free massless theory without zero mode contributions. We observe that for 
small quark masses, the contribution of the zero-mode part of the propagators completely 
saturates the propagator and is more than an order of magnitude larger than the free con- 
tribution. On the other hand, as the quark mass increases, the zero-mode contribution is 
gradually suppressed and for m q ~ 200 MeV it is only a few times larger than the free 
contribution. 

As the quark mass gets larger and larger, the quark loop contribution to QCD correlation 
functions becomes more and more suppressed. Interestingly, in the IILM, the scale m* also 
determines the regime where the quenched approximation is expected to become reliable. In 
fact, from Eq.s (fl~2"]) . (TT3T) and (1131) . we see that for quark masses m q 3> max[T/j] ~ m*, the 
contribution of the overlap matrix elements to the fermionic determinant becomes negligible. 
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FIG. 8: Sum of the zero- mode and free contributions to the pseudoscalar point-to-point correlation 
function computed in the IILM at r = 1 fm for different quark masses (diamonds) compared with 
the free contribution alone (squares). 

On the other hand, in this model, the mass contribution to the ultraviolet sector of the 
spectrum factorizes and therefore cancels out in all ensemble averages. 

The numerical results presented in the previous section reveal two interesting features 
of the Instanton Model: (i) for m q < m* the density of quasi zero-modes increases as the 
quark mass decreases and (ii) for m q > m* the near zero-mode part of the spectrum becomes 
practically independent on the quark mass. These two non-trivial dynamical effects are in 
fact related and can be explained as follows. Quark loops are known to generate strong non- 
local correlations between pseudoparticles of opposite topological charge. Such correlations 
tend to suppress configurations in which one or more pseudoparticles are located far from 
all others. Hence, the contribution of the fermionic determinant leads to a reduction of 
the density of nearly zero-modes, as completely isolated instantons are known to have exact 
zero-modes. As the quark mass gets larger, such a topological screening becomes less and 
less effective and the population of near exact zero-modes increases, explaining the rise of 
the peak of eigenvalue density near the origin. On the other hand, we have seen in section [III 
that for m q > m*, the contribution of the fermionic determinant is suppressed and fermion- 
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induced topological screening completely disappears. As a result, the low-virtuality sector 
of the Dirac spectrum stops depending on the quark mass. These are interesting predictions 
of the instanton liquid model that can be checked with unquenched lattice calculations. 

We conclude this section by emphasizing that the existence of a single mass scale m* 
at which chiral symmetry breaking begins to diminish and at which the quenched approx- 
imating begins to become valid is a characteristic prediction of the instanton model. In 
principle, in full QCD we expect the scales for such transitions to be implicit functions of 
the only dimensional parameters of the theory, Aqcd an d m q and there is no a priori reason 
for which they should coincide, unless they are driven by a common dynamical mechanism. 



VII. CONCLUSIONS 



In this work we have used the IILM to investigate the dynamical mechanisms that drive 
the dynamics and the structure of hadrons in the light quark sector of QCD. To this end, we 
have checked that the model's predictions are consistent with QCD in the pion mass regime 
< 500 MeV, where we expect chiral dynamics to play an important role. By computing the 
nucleon mass and pion mass at different values of the quark mass we have shown that the 
IILM is consistent with the existing lattice data in the range m n ~ 300 — 600 MeV. We 
have observed that the best agreement with lattice QCD data in this region is obtained for 
a value of the average instanton size p = 0.32 fm, which is close to Shuryak's earlier phe- 
nomenological estimate. On the other hand, the calculated instanton density n 3 fm -4 is 
considerably larger than in earlier phenomenological estimates and closer to lattice results. 
By studying the dependence on the quark mass of the density of eigenvalues of the Dirac 
operator, we have shown that the model also contains the correct dynamics to reproduce 
predictions of xpt in the small quark mass regime where such a theory is applicable and 
finite-order calculations are reliable. We have shown that, in the chiral limit, our simulations 
converge to the Nf = 2 prediction: lim^o = 0. We have studied the modification of 
the Dirac spectrum induced by small quark masses and found that they are proportional to 
m"A 1_a , a functional form that generates a l/m q power-law divergence in the quark mass 
dependence of the scalar three-point correlators as predicted by xpt |22J]. We have discussed 
how the qualitative structure of such mass corrections has a simple physical interpretation 
in terms of quark-loop-induced instanton-antiinstanton correlations. It would be very in- 
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teresting to check if the corrections predicted in the IILM are observed in full lattice QCD 
simulations. 

From a chiral extrapolation of the pion mass, we computed the effective chiral parameters 
f n and B and showed that they are comparable to those in QCD. Similarly, from the 
agreement of our results for the nucleon mass at different quark masses with the available 
lattice simulations, we conclude that the parameters for the baryon chiral perturbation 
theory are consistent with those extracted from QCD, although our present statistics do not 
allow us to perform a direct chiral extrapolation. 

Having checked that chiral dynamics is correctly incorporated in the instanton model, 
we have exploited our analytic understanding of instanton-induced correlations to study the 
dynamical mechanisms involved in the transition into the chiral regime at a microscopic 
level. We have identified a mass scale m* = 80 MeV above which we do not expect QCD 
correlators to be dominated by chiral dynamics and above which quenched calculations 
should begin to approximate full QCD. 

If the computational technology is developed to the point that the present statistical and 
systematic errors on the nucleon mass are significantly improved, this framework could rep- 
resent a complementary tool to perform direct chiral extrapolation of a variety of observables 
of interest, such as moments of DIS and generalized parton distributions, magnetic moments 
and form factors. Another interesting development would be to investigate if unstable light- 
quark resonances that receive contribution at order n 2 level in the instanton density, such 
as the p-meson and A-isobar, can also be well described by the 't Hooft interaction. In this 
context it would also be interesting to investigate the effect of confinement on the hadron 
masses, extending the pseudoparticle ensemble to include both regular and singular gauge 
instantons j^ . 
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